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g : ABSTRACT 

(N : 

'. The jet properties of Mrk 501 and 3C 279 are derived by fitting the broadband 

pL| . spectral energy distributions with the lepton models. The derived 7b (the break 

t:J- ! Lorenz factor of the electron distribution) are 10^ — 10^ for Mrk 501 and 200 ~ 

! 600 for 3C 279 and the magnetic field strength {B) of Mrk 501 is usually one 
order of magnitude lower than that of 3C 279, but their Doppler factors (6) are 

i-jH ! comparable. A spectral variation feature that the peak luminosity is correlated 

^ ! with the peak frequency, which is opposite to the blazar sequence, is observed 

CLi! in the two sources. We find that (1) the peak luminosities of the two bumps 

O ! in SEDs for Mrk 501 depend on 7^ in both the observer and co-moving frames, 

! but they are not correlated with B and 5; (2)the luminosity variation of 3C 

, 279 is dominated by the external Compton (EC) peak and its peak luminosity 

! is correlated with % and 6, but anti-correlated with B. These results suggest 

^ , that 7b may govern the spectral variation of Mrk 501 and 6 and B would be 

^ i responsible for the spectral variation of 3C 279. The narrow distribution of 7b 

00 ! and the correlation of 7b and B in 3C 279 would be due to the cooling by the EC 

^ i process and the strong magnetic field. Based on our brief discussion, we propose 

^ I that this spectral variation feature may be originated from the instability of the 

^ ' corona but not from the variation of the accretion rate as that for the blazar 

^ I sequence. 
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1. Introduction 

Blazars, a sub-sample of active galactic nuclei (AGNs), are composed of BL Lac objects 
(BL Lacs) and flat spectrum radio quasars (FSRQs). Their broadband spectral energy 
distributions (SEDs) are usually bimodal. The bump in the IR-optical-UV band can be 
explained with the synchrotron process of relativistic electrons accelerated in the jets, and the 
bump in the GeV-TeV gamma-ray band may be due to the inverse Compton (IC) scattering 
of the same electron population. The seed photon fields for BL Lacs are thought to be the 
synchrotron radiation themselves (the so-called SSC model, Maraschi et al. 1992; Ghiscllini 
et al. 1996; Zhang et al. 2012a), but the contributions of inverse Compton scattering by 
some external radiation fields (EC), such as the broad line region (BLR), usually dominate 
the IC peak of the SEDs for FSRQs (Dermer et al. 1992). 

Significant variability in multi-wavelengths with spectral shift is observed in blazars 
(Massaro et al. 2008; Tramacere et al. 2009; Zhang et al. 2012a) and its physical reason 
is still uncertain (e.g., Zhang et al. 2012a). Associated flux variations in the X-ray and 
gamma-ray bands were observed in some sources (Takahashi et al. 2000; Sambruna et al. 
2000). Recently, Zhang et al. (2012a) analyzed the SEDs of GeV-TeV BL Lacs and found 
that these SEDs are well modeled with the one zone Syn-|-SSC model. They showed that the 
observed luminosity variations are usually accompanied with the shift of peak frequencies of 
the SEDs (see also Tavecchio et al. 2001; Blazejowski et al. 2005) and the spectral shift may 
be related to the electron acceleration in the jet. This is consistent with that reported by 
Massaro et al. (2008), Tramacere et al. (2009), and Fossati et al. (2008), who suggested that 
the correlation between the SED peak energy of synchrotron radiation and the corresponding 
peak flux would be a signature of the acceleration or injection of the relativistic electrons at 
the higher energy end by stochastic acceleration mechanism. However, some "orphan" flares 
detected in one band without counterparts in other bands were also observed (Krawczynski 
et al. 2004; Blazejowski et al. 2005). It is uncertain that these flares are from different 
radiating regions or due to the temporal evolution of the parameters for the emitting region 
(Fossati et al. 2008). 

Well-sampled SEDs in different episodes for a given blazar are good probes for the jet 
properties since the differences in both the jets and their environments among sources should 
be eliminated. This work dedicates to investigate the spectral variations in individual blazars 
using broadband SEDs with significant variation in different episodes. We compile the SEDs 

for Mrk 501 (a typical BL Lac) and 3C 279 (a typical FSRQ). The data are presented 
in §2. They are fitted with the syn-|-SSC or syn+SSC+EC models. Our models and the 
derived parameters are described in §3. The possible relations among parameters in both 
the observed and co- moving frames and the corresponding implications are presented in §4. 
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Discussion on the spectral variation in an individual blazar and the blazar sequence is given 
in §5. A summary of our results is presented in §6. 

2. Data 

We compile the SEDs of Mrk 501 and FSRQ 3C 279 from literature. We get four well- 
sampled SEDs for both of them as described bellow. They are shown in Figure [TJ We mark 
the SEDs as "a" , "b" , "c" , and "d" states according to their peak luminosity of synchrotron 
radiation. 

Mrk 501. It is a TeV source with z = 0.034, which was identified by Whipple during 
an observation campaign of 66 hr. An average flux of (8.1 ±1.4) x 10~^^cm~^s~^ above 300 
GeV and flux variation on timescale of days was observed (Quinn et al. 1996). From 1997 
April to 1999 June, this source was detected quasi-simultaneously with BeppoSAX and the 
Cerenko Array. The observations show that the peak frequency of the synchrotron emission 
shifts from 100 keV back to 0.5 keV, and correspondingly the flux decreased (Tavecchio et 
al. 2001). The SEDs "a" and "b" of this source are taken from the broadband observations 
in 1997 April 16th and 1997 April 7th. The data of the SED "c" is taken from a 4.5 month- 
long multi-wavelength observation campaign for Mrk 501 performed from 2009 March 15 to 
August 1 (Abdo et al. 2011)0. The SED marked as "d" is taken from a multi- wavelength 
campaign with Suzaku and MAGIC performed in 2006 July. The average VHE flux above 
200 GeV in this campaign is ~ 20% of the Crab flux, with a photon index 2.8 ± 0.1 from 80 
GeV to 2 TeV (Anderhub et al. 2009). The data in the GeV band for this SED are taken 
from Fermi/L AT observations (Abdo et al. 2009). 

3C 279. The redshift of this sources is 0.536. It is the first bright 7-ray blazar reported 
by the EGRET instrument aboard the Compton Gamma-Ray Observatory (Wehrle et al. 
1998). It is also the first confirmed TeV FSRQ (Albert et al. 2008). Its SED "a" is 
obtained from a monitoring program with the Compton Gamma Ray observatory, the Ginga 
X-ray satellite, and instruments in other bands performed in 1991 June during a 7-ray flare 
(Hartman et al. 1996). The SED "b" is from observations in 2006 February (Berger et al. 
2011). The SED "c" and "d" are obtained during this source experienced a sharp gamma-ray 
flare and an isolated X-ray flare reported by Abdo et al. (2010), respectively. 



^Note that the TeV data of this SED has been corrected for the extra-galactic hght (EBL) absorption 
witli a model proposed by Franceschini et al. (2008), so wc do not make fm'ther EBL correction in our SED 
fit. 
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3. Model and SED Fits 

As shown in Zliang et al. (2012a), tlie external photon fields from the BLRs can be 
neglected in comparison with the synchrotron radiation photon fields for BL Lacs. Therefore, 
the syn+SSC model is used to fit the SEDs of Mrk 501. For FSRQ 3C 279, the contribution 
of external photon field from the BLR cannot be neglected, and thus the syn+SSC+EC 
model is used to fit the SEDs of 3C 279. The total luminosity of the BLR for 3C 279 is taken 
from Celotti et al. (1997) and the size of the BLR is calculated using the BLR luminosity 
with equation (23) presented in Liu & Bai (2006). The radiation from a BLR is taken as a 
black body spectrum and the corresponding energy density measured in the co-moving frame 
IS t/gLR = 5.37 X lO-^F^ erg cm , where the Lorentz factor of the jet is taken as its Doppler 
boosting factor, i.e., F = 5, because the relativistic jets of blazars are close to the line of 
sight. The minimum variability timescales are taken as At = 1 hr for Mrk 501 (Tavecchio 
et al. 2001) and At = 8 hr for 3C 279 (Wehrle et al. 1998). The energy distributions of 
the relativistic electrons in the radiating regions are assumed to be a broken power-law with 
indices pi and p2 as well as a break at 7b within the energy range [7min5 7max]- The details of 
the models and the strategy for parameter constraints are presented in Zhang et al. (2012a). 

The SEDs of Mrk 501 and 3C 279 are well fitted with the models. Our best SED 
fits and the contours of the probability p {p oc e'^l"^) for the parameters B and 5 at 1 cr 
significance level are also presented in Figure [U The contours are thin and slanted ellipses, 
especially for Mrk 501, which are consistent with that reported by Zhang et al. (2012a). 
The EC components of SEDs for 3C 279 present further constraint on 5, hence make tighter 
constraint on B and 5 than that for Mrk 501. The parameters of the best fit to the SED as 
shown in Figure 1 are marked with stars in the contours. We also calculate the bolometric 
luminosities (Lboi) with our best SED fits. Our results are reported in Table 1. One can 
observe that 7b ranges in 10^ - 10^ for Mrk 501, but it is 200 ~ 600 for 3C 279. The values of 
5 for the two sources are comparable within errors. The magnetic field strength of Mrk 501 
is roughly one order of magnitude lower than that of 3C 279. 

4. Observed Spectral Variation and the Possible Physical Origins 

4.1. Spectral Variation of Synchrotron Radiation Bump 

The observed peak luminosity [L^ of the synchrotron radiation as a function of the 
corresponding peak frequency (z/g) for the SEDs in our sample is shown in Figure [2] (a). By 
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making the beaming factor correctioiiG for and Lg, we show the relation of Ug — in the 
CO- moving frame in Figure |2] (b), where the prime denotes the data in the co-moving frame. 
To study the dependence of Lg and L'^ on the physical parameters of the jets, we show 
(or Lg) as functions of B, 6, and 7b in Figure [3l The results of pair correlation analysis with 
the Pearson method are reported in Table 

• Mrk 501. A clear correlation between Lg and Ug in the observer frame is found, i.e., Lg oc 
^o.45±o.i5^ which is consistent with that observed in Mrk 421 {Lg oc z/0-42±o.o6^ Tramacere 
et al. 2009). L'^ is also correlated with z/g, but the slope changes to 0.85 ± 0.23. No 
clear dependence of Lg and L'^ on B and 5 is observed. However, both Lg and L'^ are 
correlated with 7b, i.e., Lg oc ^o-87±o.30 ^ ^i.7o±o.7i^ These results suggest that 
the observed Ug — Lg correlation could be dominated by the variation of 7b. 

• 3C 279. The variations of Lg and L'^ of 3C 279 are much smaller than that of Mrk 501. 
Although a tentative correlation between Lg and Ug is also observed, the slope of the 
correlation is much shallower than that in Mrk 501, i.e., Lg oc i/0-28±0-i9. Even more, we 
do not find any correlation between Lg and z/g. This is reasonable since 7b of 3C 279 
in different states is in a very narrow range. We do also not find any dependence of Lg 
on 7b, B, and 6 as shown in Figure [31 Tentative 7b — Lg correlation and B — Lg anti- 
correlation are found in 3C 279. We infer that they may be from an anti-correlation 
between B and 7b, as shown in Figure HI A clear anti-correlation of S — 7b is observed 
in 3C 279, which is 7b oc i^-i-66±o.i9 rjj-^g observed anti-correlation of i? — 7b in 3C 
279 is consistent with the radiation cooling for the electrons by the magnetic field. 

Similar z/g — Lg correlation was also found by Massaro et al.(2008) for some high-frequency- 
peaked BL Lacs (see also Zhang et al. 2012a). They proposed that this correlation is a 
signature of synchrotron radiation. Under the mono-frequency approximation condition, we 
get Lg oc N'y^B'^6'^ and z/g oc 'y^B6, where N is the number of emitting electrons. The 
comparison between the observed correlations and the model predictions can provide some 
information about which parameter may dominate the observed correlations. One can infer 
that L' oc z/' if 7b dominates the relation of z/' — L' and L' oc u'!^ if B dominates the relation 
of z/g - Lg (see also Massaro et al. 2008). We find that Lg oc i/^o.85±o.23 y^^y^ ^q]^^ indicating 
that the relation of z/g — Lg in Mrk 501 should be due to the variation of 7b. However, no 



^Here we do not consider the errors of S and just use the central value of S to correct the parameters. 

■^Considering the small sample statistics and large errors of the parameters, we also use a bootstrap 
method to examine the correlation coefficients of these correlations. The results are consistent with that of 
the Pearson correlation analysis and are also shown in Table 2. 
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significant spectral variation of synclirotron process bump is observed for 3C 279, especially 
in the co-moving frame. 

4.2. Spectral Variation of the IC Bump 

We present similar analysis for the IC bump as that for the synchrotron radiation 
bump. The IC bump is dominated by the SSC process for Mrk 501 and the EC process 
for 3C 279, therefore, we give the analysis of the SSC and EC bumps for Mrk 501 and 3C 
279, respectively. Figure shows the correlations of — in the observed and co-moving 
frames and Figure |6] displays Lc (or L'J as functions of B, 6, and 7b for the two sources. The 
Pearson correlation analysis results are also reported in Table 2. 

• Mrk 501. It is found that Lc is correlated with Uc in both the observer and co- moving 
frames. The best fits give Lc oc i/0-85±o.i4 ^^^^ ^ ^^i.46±o.o3_ rjj-^g power-law index of 
the observed i>c — Lc relation is much larger than that of the — Lg relation, especially 
in the co-moving frame. This is inconsistent with the SSC model. Note that the IC 
peak at the TeV gamma-ray band would be affected by the Klein-Nishina (KN) effect, 
which would result in a larger slope in the observed Uc — Lc relation than that in the 
z/g — Ls relation, especially for a transition from the Thomson regime to the KN regime 
(Tramacere et al. 2011). The IC peak of SED "a" for Mrk 501 is in the Thomson 
regime and the other three SEDs are in the KN regime. Therefore, the discrepancy of 
the slopes in the IC process and the synchrotron radiation bump may be due to the 
KN effect. Both Lc and L'^ are correlated with 7b. The best fits give Lc oc ^o-86±o.20 
and L'^ oc ^i-69±o-59_ -^^ clear dependence of Lc and on B and 6 is found. Therefore, 
the variation of the IC peak for this source should be also due to the variation of 7b 
as that for the synchrotron peak. 

• 3C 279. Its IC peak is dominated by the EC process. Its Uc is also correlated with 
Lc in both the observer and co-moving frames. The best fits give Lc oc z/i o6±o.22 

L'^ oc zy^o-55±o.23_ j^^^g ^^Y^aX the slopes of relation Vc — Lc in 3C 279 and Mrk 501 are 
roughly consistent within the error bars, but the slope of v'c — L'^ in 3C 279 is much 
smaller than that in Mrk 501. This is because the IC bump of 3C 279 is dominated by 
EC process and the select SEDs of 3C 279 at gamma-ray band are not affected by the 
KN effect. Lc is strongly dependent on 5, i.e., Lc oc (56.2o±i.70^ jg ^^^^^ indicating 

that the Doppler boosting factor is a great ingredient for the variation of Lc- Both 
Lc and are also anti-corrected with B, i.e., Lc oc i^-4-42±o.46 ^ ^-i.96±o.68_ 

This implies that the magnetic field may also play an important role in the spectral 
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variation of 3C 279. Similar to Mrk 501, both and are corrected with 7^ for 3C 
279, i.e., Lc oc -)/2-55±o.54 ^'^ ^ ^i.i9±o.37_ rpj^jg suggests that 7b is also a key factor 

for the variation of the EC peak, and is anticorrelated with B, being consistent with 
the radiation cooling for the electrons by the magnetic field. 

As shown above, the variation of the IC peak for Mrk 501 is also dominated by the electron 
spectrum, similar to its synchrotron peak. The variation of the IC peak for 3C 279 depends 
on 5, B, and 7b. The electron spectrum is critical since both the synchrotron radiation and 
the IC process are highly dependent on it. Significant variabilities in both the synchrotron 
process and IC bump of Mrk 501 are related to the variation of 7b, which ranges within 
10^ ~ 10^, implying that the particle acceleration in different epoches of Mrk 501 would be 
different. Note that the environments of the emitting regions of BL Lacs and FSRQs are 
dramatically different. The radiating region of 3C 279 is suggested to be inside the BLR. The 
electrons thus would be cooled down by the EC process. On the other hand, the magnetic 
field strength in 3C 279 is usually one order magnitude higher than that in Mrk 501 and the 
electrons may also be cooled through the radiation cooling effect by the magnetic field. These 
effects would result in a narrow distribution of 7b (i.e., 7b = 200 ~ 600) among different 
states for 3C 279, and thus there is no significant variation in the synchrotron bump. The 
variation of the IC peak for 3C 279 is dominated by the EC process. The energy density 
of the external photon field is magnified by (~ 5^) and the energy of seed photon is 
magnified by F (~ 5) in the co-moving frame, thus a small variation of 5 would result in 
significant variations of L[. and Hence, the variation of IC peak for 3C 279 is dominated 
by the variation of 5. However, the variation of B may be also responsible for the variation 
of IC peak, which result in a correlation of 7b — (or 7b — L'^) because of the cooling effect. 



5. Discussion 

It was found that high-luminosity FSRQs tend to have a low peak frequency and low- 
luminosity BL Lacs tend to have a high peak frequency. This is the so-called blazar sequence, 
which may be related to the different environments of emitting regions for different types 
of blazars (e.g., Ghisellini et al. 1998). With the discovery of "blue quasars", Ghisellini & 
Tavecchio (2008) further proposed that, more physically, the blazar sequence is due to the 
difference of the black hole (BH) masses and accretion rates among these sources. In this 
scenario, the "blue quasars" should have a large BH mass and an intermediate accretion 
rate. Their emission regions are beyond the BLR. Moreover, the red low-luminosity blazars 
should exist, which may have a small BH mass and a relatively large accretion rate. Chen 
& Bai (2011) reported that the properties of the narrow fine Seyfert 1 galaxies, which have 
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a low-peak-frequency and low- luminosity, are similar to this kind of blazar. As presented in 
Figures [2] and |5l the data points of 3C 279 are in the left-top and the data points of Mrk 
501 are in the right-bottom of the figures, illustrating the feature of the blazar sequence. 
The observed z/g — correlation in individual blazars is opposite to the blazar sequence and 
is referred as blazar anti-sequence of spectral variation in individual blazars (Zhang et al. 
2012b). 

The results in this analysis indicate that the environment effect may be also fundamental 
for the observed different spectral variations of 3C 279 and Mrk 501. The significant differ- 
ences of the electron energy distributions between Mrk 501 and 3C 279 may be mainly due 
to the cooling effect of the external photons from the BLR of 3C 279. However, the spectral 
variation among different states of an individual source is independent of its environment. 
As shown in our analysis the variation of the electron spectrum may result in the spectral 
variation of BL Lacs, and the bulk motion of the emitting region as well as the magnetic 
field strength would be responsible for the spectral variation of FSRQs. 

Different from the blazar sequence, the rapid spectral variation in individual blazars 
would be not originated from the variation of the accretion rate. The variability timescale 
produced by the variations of the accretion rate can be estimated by the viscous timescale 
of disk, 

~ 3 X 10'a-'/'M^,'^''my'R%' s, (1) 

where Rw = R/{W'^ cm), Miq = M/(10^^ g s"^), mi = M/Mq. Taking the typical values 
of those parameters for blazars, Mbh = 10^'^Mq, M = O.lMEdd, R = ^ORs, and a = 0.1, 
we get tvisc = 6.03 x 10® days, which is much larger than the observed variability timescales 
of blazars. The instability of the accretion disk cannot produce the rapid variability as that 
observed in blazars. On the other hand, the variation of the accretion rate may result in 
the luminosity variations of both the disk and the broad emission lines. We cannot test this 
scenario with Mrk 501 since BL Lacs have no or very weak emission lines and their disk 
luminosity is also very weak, which may be much lower than the non-thermal radiation from 
the jet. It is found that the strength of Lya of 3C 279 is nearly constant when the continuum 
flux varies, as that also observed in 3C 273 (Koratkar et al. 1998). Therefore, the blazar 
anti-sequence of spectral variation in individual blazars would be not originated from the 
change of the accretion rate. 

The rapid variability of blazars should be dominated by the jet emission. It is generally 
believed that jet launching is connected with the accretion disk and corona in a source 
(Armitage & Natarajan 1999; Merloni & Fabian 2002; Cao 2004, Chen et al. 2012). Using 
the observations of BeppoSAX, Grandi & Palumbo (2004) untangled the emissions of jet, 
accretion disk, and corona for 3C 273, and found that the emission of corona becomes weak 
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when the jet emission increases in X-ray band, implying that the variation of the jet emission 
is related with the variation of the corona emission. Therefore, the instability of the corona 
may result in the variation of the physical condition of jets and lead to the variation of the 
jet emission. On the other hand, in a BH X-ray binary, the X-ray spectrum becomes soft 
when the corona-dominated emission drops, and the properties of its jet, such as density and 
magnetic field, is also changed (Fender & Belloni 2012). This is consistent with our results, 
implying that the jet emission of blazars may be also related with the corona emission. The 
dominated mechanism of corona variation may be magnetic reconnections, like the corona 
activity in the Sun (Zhang et al. 2000; Zhang 2007). For a super-massive black hole, the 
timescale of magnetic reconnections in the disk is from hours to days, which is consistent 
with the observed variability timescale of blazars. 



6. Summary 

SEDs observed at four epochs for Mrk 501 and 3C 279 are compiled from literature to 
investigate the spectral variation in individual blazars. Our results are summarized below. 

• The SEDs of Mrk 501 and 3C 279 are well fitted with the syn+SSC and syn+SSC+EC 
models, respectively. The derived values of 75 are 10^ — 10^ for Mrk 501 and 200 ~ 600 
for 3C 279. The magnetic field strength of Mrk 501 is usually one order of magnitude 
lower than that of 3C 279. The Doppler factors of the two sources are comparable. 

• For Mkn 501, the peak luminosities of the synchrotron radiation and IC bumps are 
tentatively correlated with the corresponding peak frequencies and they are also corre- 
lated with 76 in both the observer and co-moving frames, but they are not correlated 
with B and S. 

• For 3C 279, no significant variation is observed for the synchrotron radiation peak and 
its luminosity variation is dominated by the EC peak. Its EC peak luminosities in 
both the observer and co-moving frames are correlated with the peak frequencies, are 
also correlated with 7b and S, but anti-correlated with B. 

These results suggest that the variation of the electron spectrum may result in the 
spectral variation of Mrk 501, and more complicated, the variations of the electron spectrum, 
the bulk motion of emitting region, as well as the magnetic field strength would be responsible 
for the spectral variation of 3C 279. However, the cooling effect by external field photons may 
smear out the difference of electron spectrum in different states of 3C 279 and its spectral 
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variation would be dominated by the EC component, hence the magnification of the external 
photon field by the bulk motion of the jet is an essential reason for the spectral variation. 

Our results show that the spectral variation in individual blazars is different from the 
blazar sequence for different types of blazars. We refer this spectral variation feature as blazar 
anti-sequence of spectral variation in individual blazars. Based on our brief discussion, we 
propose that this spectral variation feature may be originated from the instability of the 
corona but not from the variation of the accretion rate as that for the blazar sequence. 
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Table 1. Observations and SED fit results for Mrk 501 and 3C 279 



Source 


State^ 


Pi 


P2 




7b 


'Ymax 




6 


B 

(G) 




Lhoi 

(crg/s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Mrk 501 


a 


1.86 


4.6 


2 


9.65E5 


1.5E8 


1E4 


1529 


0.4^:^3 


10.8 


2E46 




b 


2 


3.2 


20 


2.56E5 


9E6 


2.2E4 






9.9 


5.5E45 




c 


2.24 


3.72 


300 


3.56E4 


1E6 


3.9E4 


oc47 
^Dl6.5 




11.6 


1.7E45 




d 


2.64 


3.3 


900 


1.16E5 


8E5 


1.7E6 


0941 


n 1 30-42 

'-'•^^"'0.07 


11.1 


1.9E45 


3C 279 


a 


2.2 


4.2 


2 


451 


3E3 


1.3E5 


2i?i 


3.4p 


5.5 


4.7E48 




b 


2.2 


4.2 


2 


563 


7E4 


1.2E5 


2o.5fi:i 


3-2ii 


13.1 


3.7E48 




c 


2.2 


4.46 


2 


394 


3E3 


5.8E4 


9922.8 

''^20.8 


3.654;! 


4.6 


2.3E48 




d 


2.1 


4.3 


45 


223 


1E4 


2.1E5 


i5.5}i:i 


5.4i3 


4.3 


5.6E47 



'^The state of each SED, which is defined according to the peak luminosity of the synchrotron 
radiation bump. 

'^It is calculated using equation (23) given in Tavecchio et al. (1998). 

Note. — Columns: (3) (4) The energy indices of electrons below and above the break; (5) (6) (7) 
The minimum, break, and maximum Lorentz factors of electrons; (8) The electron density parameter 

Nq; (9) The beaming factor S; (10) The magnetic field strength B: (11) The lower limit of ^ derived 
by the 7-ray transparency condition; (12) The bolometric luminosity of each SED. 
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S 6 

Fig. 1. — The observed SEDs [scattered data points) with model fitting [lines) for Mrk 501 
and 3C 279, and the contours of 1-a significance level in the B-6 plane. The four SEDs are 
named as states "a" [blue points and thick solid line), "b" [magenta points and dashed line), 
"c" [red points and thin solid line), and "d" [green points and dotted line), respectively. The 
red stars stand for the best fit parameter sets of 6 and B. 
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Table 2. Results of the best linear fits with the Pearson correlation analysis method for the 
parameter sets (x — y) of Mrk 501 and 3C 279. The fitting function is logy = a + 61ogx. The 
Spearman correlation coefficient is r and the chance probability is p. 



Mrk 501 3C 279 







a 


b 


r 


P 


^ bootstrap 


a 


b 


r 


P 


^ bootstrap 




-Ls 


36.6±2.8 


0.45±0.15 


0.90 


0.10 




43.0±2.6 


0.28±0.19 


0.72 


0.23 




I's 




25.3±3.8 


0.85±0.23 


0.94 


0.06 














7b 


-Ls 


40.1±1.6 


0.87±0.30 


0.90 


0.10 




45.1±0.7 


0.69±0.27 


0.87 


0.13 




7b 


-l's 


30.4±3.8 


1.70±0.71 


0.86 


0.14 


















22.1±3.6 


0.85±0.14 


0.97 


0.03 


0.90 


22.5±5.2 


1.06±0.22 


0.96 


0.04 


0.96 




-< 


2.8±0.9 


1.46±0.03 


~ 1.0 


5.6E-4 


0.96 


30.1±5.1 


0.55±0.23 


0.86 


0.14 


0.83 


B 














50.0±0.3 


-4.42±0.46 


-0.99 


0.01 


-0.87 


B 


-< 












43.4±0.4 


-1.96±0.68 


-0.90 


0.10 


-0.67 


S- 














39.4±2.2 


6.20±1.70 


0.93 


0.07 


0.98 


7b 


-ic 


39.7±1.0 


0.86±0.20 


0.95 


0.05 


0.95 


40.8±1.4 


2.55±0.54 


0.96 


0.04 


0.96 


7b 


-L'c 


30.0±3.1 


1.69±0.59 


0.90 


0.10 


0.90 


39.2±1.0 


1.19±0.37 


0.92 


0.08 


0.88 


B 


- 7b 












3.6±0.1 


-1.66±0.19 


-0.99 


0.01 


-0.91 



^The correlation coefficient rbootstrap is obtained by the bootstrap method as described in the text. 
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Fig. 2. — Peak frequency of synchrotron radiation as a function of the peak luminosity in 
the observer {Panel a) and co- moving {Panel b) frames. The best fit lines and the slopes are 
also shown in each panel. 
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Fig. 3. — Magnetic filed strength B, the beaming factor 6, and the break Lorentz factor of 
electron 7b as a function of the peak luminosity of synchrotron radiation in both the observer 
(Ls) and co-moving (Lg) frames for Mrk 501 and 3C 279. 
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Fig. 4. — Break Lorentz factor of electron 7b as a function of the magnetic filed strength B 
for Mrk 501 and 3C 279. 




Fig. 5. — Peak frequency of IC process as a function of the peak luminosity in the observer 
[Panel a) and in the co-moving [Panel b) frames. 
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Fig. 6. — Magnetic filed strength B, the beaming factor 6, and the break Lorentz factor of 
electron 7b as a function of the peak luminosity of the IC bump in both the observer (Lc) 
and co-moving (LJ frames for Mrk 501 and 3C 279. 



